M1O6  830 


INVESTIGATION  OF  GROUND  SHOCK  EFFECTS 
IN  NONLINEAR  HYSTERET1C  MEDIA 


Report  4 


EFFECT  OF  A  STEP  LOAD  MOVING  WITH  CONSTANT 
SUPERSEISMIC  VELOCITY  ON  A  HALF-SPACE  OF  A 
VARIABLE  MODULUS  MATERIAL 


by 

A,  T.  Matthews,  I.  Sandler,  H.  H.  Bleich 


March  1970 

Sponsored  by  Defense  Atomic  Support  Agency 

Conducted  for  U.  S.  Army  Engineer  Waterways  Experiment  Station,  Vicksburg,  Mississippi 
Under  Contract  No.  DACA39-67-C-004Q 

By  Paul  Weidlinger,  Consulting  Engineer,  New  York,  New  York 


ARMY'MflC  VICKSBURG.  HISS. 


This  document  has  been  approved  for  public  release  and  sale;  its  distribution  is  unlimited 


t 


FOREWORD 

The  current  work  on  variable  moduli  models  is  part  of 
Contract  No.  DACA39-67-C-0048 ,  "Investigation  of  Ground  Shock 
Effects  in  Nonlinear  Hysteretic  Media",  being  conducted  for 
the  U.§.  Army  Engineer  Waterways  Experiment  Station  (WES)  under 
DAS A  sponsorship. 

The  contract  is  being  monitored  by  Dr.  J.S.  Zelasko, 
Impulse  Loads  Section,  Soil  Dynamics  Branch,  under  the  general 
supervision  of  Mr.  J.P,  Sale,  Chief,  Soils  Division,  WE£. 

COL  Levi  A.  Brown,  CE,  is  the  contracting  officer. 


ABSTRACT 


The  report  considers  the  two  dimensional  effects  of 
a  step  wave  progressing  with  constant  superseismic  velocity 
on  the  surface  of  a  half-space.  The  material  treated  is 
isotropic  and  it  is  assumed  that  incremental  relations  between 
both  deviatoric  and  volumetric  stress  and  strain  depend  not 
only  on  the  instantaneous  values  of  those  quantities,  but 
also  on  bulk  and  shear  moduli  which  differ  according  to 
whether  initial  loading,  unloading,  or  reloading  occurs. 

Simple  closed  form  solutions  are  obtained  when  the  moduli 
K  and  G  are  constant  and  lie  within  certain  limits.  For  the 
more  general  case,  when  K  and  G  are  functions  of  the  first  and 
second  invariants  of  stress,  solutions  requiring  only  quadra¬ 
tures  are  found. 

As  an  auxiliary  study,  a  problem  involving  a  naif-space 
of  fluid  with  a  bilinear  pressure  volume  relation  is  solved. 
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LIST  OF  SYMBOLS 


u,  v 

V 

V  ,  V 

x  ’  y 
X 

x  >  y «  2  s  ?  >  ri  i  £ 

V 

Aa,  Av,  etc. 


ekk 
0,  5 


Velocity  of  propagation  of  elastic 
P-waves  and  S-waves. 


'ij 


Strain  deviators. 

Shear  and  bulk  moduli,  respectively. 

Constants  used  in  express  ions  for  bulk 
and  shear  moduli. 

Invariants,  Eqs .  (5)  and  (6). 

Nondi mens ional  pressure. 

Stress  deviators  with  respect  to  axes 
x,  y,  etc. 


Time  . 

N'ondimensional  velocities. 

Velocity  of  surface  load. 

Velocities  with  respect  to  axes  x  and  y. 
Nondimensional  quantity  defined  by  Eq .  (17), 
Coordinates  defined  in  the  text. 

Shear  strain. 

Increments  of  <7,  v,  etc-,  at  a  discon- 
tinui ty . 

Volumetric  strain. 

Angles  related  to  directions  of  principal 
stress  defined  in  the  text. 


A  Per.sity. 

0  ,  o'  ,  a  ^  ,  o 2  ,  3  Stress  components  with  respect  to  axes 

y  ^  x ,  y  ,  e  t  c . 

X  Snear  stress. 

<p,  <p  ,  <p,.  ,  Position  angles  measured  from  the  surface 

1  '  an-’,  defined  in  the  text. 


*) 


Other  symbols  in  the  text  are  defined  as  they  occur. 
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INTRODUCTION. 


In  References  (1]  and  [2]  a  material  model  was  described 
which  can  be  fitted  reasonably  well  to  the  behavior  of  soils 
in  uniaxial  strain  and  triaxial  shear  tests.  The  model  was 
introduced  for  the  investigation  of  wave  propagation  and  is 
characterized  by  the  fact  that  the  stress-strain  diagrams  for 
typical  loading-unloading  cycles  show  hysteresis  loops.  This 
is  achieved  by  prescribing  different  bulk,  and  shear  moduli  in 
loading  and  unloading.  The  model  is  thus  an  alternative  to 
elastic-plastic  ones.  The  new  model  does  not  require  any 
statement  equivalent  to  a  yield  condition.  The  following  will 
only  give  details  required  in  the  present  investigation.  For 
full  information  the  reader  is  referred  to  the  above  references. 

In  order  to  gain  an  understanding  concerning  the  propa¬ 
gation  of  dynamic  disturbances  a  previous  study,  Ref.  [3], 
presented  closed  form  solutions  for  one  dimensional  propa¬ 
gation  of  plane  waves.  The  present  study  is  concerned  with 
the  plane  problem  of  the  two  dimensional  effects  of  a  step 
wave  progressing  with  supersonic  velocity  on  the  surface  of 
a  half-space.  Solutions  of  such  two  dimensional  problems 
by  ether  than  purely  numerical,  finite  difference  methods 
not  only  add  to  the  understanding  of  the  problem,  but  permit 
also  a  check  on  the  correctness  anu  effectiveness  of  numerical 
codes  for  the  multidimensional  situations.  The  desire  for 
such  a  check  was  che  prime  motivation  for  the  present  in¬ 


vest  ig  :aon. 
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The  material  treated  Is  isotropic  and  the  model  assumes 
that  the  incremental  relations  between  volumetric  stress  and 
strain,  and  between  deviatoric  stress  and  strain  depend  not 
only  on  the  instantaneous  values  of  these  quantities,  but 
also  on  the  question  of  whether  the  change  in  these  quantities 
occurs  during  initial  loading,  unloading,  or  reloading.  The 
stress-strain  relations  appear  thus  in  form  exactly  as  in  a 
conventional  elastic  material 


and 


(1) 


'ij 


2G  Si1 


i  t  j 


(2) 


where  ,  ekk  represent  the  first  invariants  of  stress 
and  strain,  respectively,  and.  ,  e^  are  the  deviatoric 

stress  and  strain  components.  The  moduli  are  expressed  in 
terms  of  the  invariants  of  the  stress  tensor.  K  is  a 
function  of  ,  while  G  may  be  a  function  of  Jj  and  J2 

K  -  K  +  k(J, )  O) 

o  1 

6  -  Go  +  g(J2  ,  Jx)  <*> 

where 

*>  0^  +  a2  +  a3  (5) 

J2  "  1  6ij  sij 


(6) 


However,  the  dependency  on  and  Jj  is  expressed  by  different 
functions  for  initial  loading,  unloading,  or  reloading. 


Consider  the  relation  between  and  where  the  value  3K 
represents  the  slope  of  the  diagram.  For  a  material 

which  initially  hardens  such  a  diagram  will  be  in  general  of  the 
form  of  Fig.  (la).  The  path  1-2  represents  initial  loading, 
while  2-3  represents  unloading  \  where  energy  considerations 
require  that  point  3  be  to  the  right  of  the  initial  path  1-2. 

On  reloading  to  point  2,  according  to  the  model,  the  material 
will  retrace  the  path  3-2  until  point  2  is  reached.  On  further 
loading,  say  point  4,  it  follows  the  initial  loading  law.  The 
simplest  case  occurs  when  the  functions  K  are  constants, 


different  for  loading  and  unloading.  The  resulting 
diagram  is  shown  in  Fig.  (lb). 
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Consider,  for  simplicity  the  situation  when  as  well  as 
the  deviators  s^  in  a  Cartesian  coordinate  system  £,r),C  are 
Kept  constant  when  l=j.  Only  one  of  the  shears  s^  =  T  is 
varied,  while  the  other  two,  s ^  =  s^  ,  vanish.  To  stud/  the 
relation  between  T  and  the  corresponding  strain  deviator  e^  '( , 

requires  simultaneous  consideration  of  the  diagrams  of  and  T , 
Figs.  (2a, b).  Initial  loading  is  represented  by  1-2,  On  un¬ 
loading,  2-3,  the  invariant  d..  will  reach  a  minimum  value  F.-r 
T j  »  0.  Further  unloading  in  T  represents  reloading,  as 
increases  again,  until  at  point  4  the  value  equai.*)  ch« 
previous  w.>:  ;.uum  J ^  .  Further  decrease  in  T,  path  4-1,  brings 


iI'imusi'  of  the  intended  application  to  soils,  the  model 
implies  that  remains  compressive. 
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the  initial  loading  law  for  G  into  force,  as  shown  in 
Figs.  (2a, b).  If  the  functions  G  are  constants,  the  x-Y 
diagram  consists  of  straight  lines,  as  shown  in  Fig.  (2c). 

It  will  be  subsequently  shown  that  simple,  nearly 
trivial  closed  form  solutions  for  the  problem  to  be  studied 
can  be  obtained  if  the  functions  K  and  G  are  constants,  and 
lie  between  certain  limits.  For  the  more  general  case  where 
K  and  G  are  functions  of  J  ^  and  solutions  requiring  only 
quadratures  will  be  found. 

Prior  to  studying  a  solid  half-space  of  the  new  model 
material  under  progressing  loads,  it  seemed  educational  to 
consider  a  similar  model  for  fluids,  with  pressure  volume 
relations  according  to  Fig.  (1).  Fcr  progressing  step  loads 
only  uninteresting  results  are  obtained,  because  in  these  cases 
the  pressure  field  becomes  simply  a  leading  shock  front  (where 
the  pressure  jumps)  followed  by  a  uniform  pressure.  The 
result  is  thus  entirely  unaffected  by  the  unloading  law.  The 
latter  would  affect  the  results,  however,  if  the  pressure 
decays,  say  linearly,  after  the  initial  sudden  rise.  Results 
in  closed  form  for  bilinear  fluids  were  obtained  and  are 
presented  in  Appendix  B.  They  may  be  a  useful  guide  in  further 
studies  concerning  solids  subjected  to  decaying  loads. 


II  FORMULATION  OF  THE  BASIC  EQUATIONS. 


Figure  (3)  indicates  the  half-space  with  a  system  of 
stationary  Cartesian  coordinates.  The  x-axis  is  in  the 
direction  of  motion  of  the  stop  load,  the  y-  and  z-axes  are 
normal  to  the  x-axis  in  and  out  of  the  plane  of  the  figure, 
respectively.  The  analysis  considers  the  case  of  plane  strain, 
Ezz  “  0,  when  the  velocity  V  of  the  step  load  is  larger  than 
the  largest  possible  characteristic  or  shock  velocity  which  is 
within  the  domain  of  the  solution. 


The  governing  equations  for  the  problem  are  in  appearance 
the  same  as  those  for  the  steady-state  two  dimensional  problem 
for  a  conventional  clastic  material.  The  equations  of  motion 
ere 


fs  .  i!^i 

3x  3  3x 


3s  ,  3J. 

— x.  +  A _ L  +  „  D 

3y  3  3y  3x  ^ 


3t 

3y 

3t 


3v 

_ > 

3t 

3v 

3 1 


(7) 

(8) 


The  constitutive  equations  are 


.  2G  *  ^Vy 

*x  +  9f  J1  =  2G  3^ 

Or  •  <*V 

*,  +  H  «  ST 


3v  3v 

^  G(3T  +  ^ 


Sx  +  ®y  "  fl  jl  "  0 


(9) 

(10) 
(11) 
(12) 


In  the  present  problem  the  moduli  G  and  K  are  functions  of  the 


invariants  as  described  in  the  Introduction. 
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Previous  experience  with  steady-state  two  dimensional 
problems  in  inelastic  materials,  Refs.  [4],  [5],  suggests  a 
procedure  for  the  solution.  In  brief,  the  follov/ing  steps 
are  required: 


1. 


L. 


3. 


4. 


Introduction  of  the  fact  that  in  a  steady-state 
problem  a'.l  unknowns  can  only  be  functions  of  the 
variable  £  -  x  -  Vt. 

Elimination  of  the  velocities  v  and  v  . 

x  y 

Dimensional  considerations,  Ref.  [4],  suggest 
the  possible  existence  of  solutions  in  which  the 
stresses  are  solely  functions  of  the  variable 
•«  cot  ^(S/y).  Its  introduction  reduces  the 
partial  'differential  equations  to  a  set  of 
ordinary  ones. 

Finally,  the  unknowns  s  ,  s  and  T  are  replaced 

x  y 

by  the  principal  stress  deviators  s^  ,  s2  and 
the  angle  0  between  the  direction  of  s^  and  the 
horizontal,  Fig.  (3), 


s  * 

X 

2 

8^^  cos  0  +  s2 

sin20 

(13) 

6  « 
y 

2 

s^  sin  0  +  s2 

2 

cos  0 

(14) 

T  «* 

(sL  -  s 2 )  sin 

8  cos  0 

(15) 

This  leads  to  the  four  differential  equations 


— 

-1 

-1 

2G 

0 

» 

8  .. 

3K 

si 

sin2\|/ 

2 

cos  ip 

G 

i-ax  -f 

-sin  2\p 

■ 

s2 

sin  2ip 

sin  2ip 

2  sin  2ip 

G 

-2  (1-2X  —) 

,  G 

sin  ip  -X  -jr 

G  , 

X  —  -  cos  <p 

-cos  2ip 

0 

( s 1  -  s2>e' 

_  _ 

(16) 


t 

where  s.  =  tt-  «  etc.,  and 
1  d<p 


X  *»  i  (£— )  sin2<(> 
1  Go 


¥ 
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(17) 

(18) 


Equations  (16)  are  linear  and  homogeneous ,  so  that  the 

i  i  i 

derivatives  of  the  unknowns  s1  ,  s 2  and  ,  and  the  value 
(Si  -  vanish  unless  the  coefficient  matrix  is  singular, 

requiring 


K  +  -  G 

(2X  -  •—)  ( 2X - - — ) 

o  o 


(19) 


In  any  region  of  finite  extent  in  ‘p ,  changes  in  stress 
or  direction  6  can  thus  only  occur  if  one  of  the  two  factors 
of  this  equation  vanishes.  Substitution  of  th  respective 
value  of  X  into  Eqs.  (16)  defines  ratios  between  the  rates 
of  change  of  the  stresses  and  that  of  6*. 


The  values  of  X  are 


* 


"  8  ~r“ 

4 

k  +  —  r. 


or 

2X  =f 

o 

When  Che  last  equation  applies,  the  invariant  and  thus  the 
value  of  the  mean  pressure  do  not  change  at  all. 


Regions  in  <j>  where  changes  of  stress  occur  can  exist 
only  if  K  +  G  and  G  are,  respectively,  functions  of  J ^ 

and/or  d ^  .  If  one  or  two  of  these  quantities  do  not  depend 
on  the  stresses,  the  above  equations  give,  instead,  one  or 
two  values  X  which,  through  Eq.  (17)  define  locations  <f)  of 
shock  fronts.  Equation  (20)  defines  the  location  of  a  dis¬ 
continuity  in  the  normal  stress  E<1»  (2D 

defines  a  location  of  a  discontinuity  in  the  shear  stress 
in  the  plane  of  the  shock  front.  Both  types  of  shock  fronts 
propagate  with  their  respective  characteristic  velocity. 

If  K  +  y  G  and  G  are  functions  of  the  stresses,  dis¬ 
continuities  propagating  with  n^ncharactaristic  velocities 
may  also  occur.  The  velocity  and  corresponding  locations  4* 
of  such  shocks  are  not  given  by  Eqs.  (20,21),  nor  defined  by 
the  differential  equations,  Eqs.  (16),  which  apply  only  in  a 
region  of  a  continuous  stress  field.  The  appropriate 
equations  must  be  obtained  for  the  nonlinear  stress-strain 
relations  in  combination  with  momentum  considerations 
(Rankine-Hugoniot  relations).  Tha  situation  is  discussed  in 
Appendix  A. 


(20) 


(21) 


Ill  BILINEAR  MATERIAL. 


In  order  to  obtain  insight  into  the  nature  of  the  stress 
field  caused  by  a  progressing  step  load,  Fig.  (3),  the  special 
case  in  which  the  functions  k  and  g  in  Eqs .  (3,4)  vanish  is 

<  insidered  first.  In  this  case  the  expressions  for  K  and  C  do 
not  depend  upon  the  stresses  but  on  whether  or  not  the  material 
is  being  loaded,  unloaded  or  reloaded.  Due  to  the  linear 
nature  of  the  stress-strain  relation,  Figs,  (lb, 2c),  this  case 
will  be  referred  to  as  "bilinear",  and  the  moduli  for  -oading 
will  be  designated  by  Kq  and  Gq  while  those  for  unloading¬ 
reloading  are  Kq  and  Gq  . 


As  mentioned  in  Section  II,  Eqs.  (20, 21)  when  applied  to 
the  bilinear  material  do  not  permit  regions  of  finite  extent 
in  <J>  within  which  the  stresses  change.  Instead,  these 
relations  define  locations  of  discontinuities  in  stress  and 
velocity,  similar  to  dilatational  (P)  and  shear  (S)  fronts  in 
a  conventional  linear  elastic  material.  The  locations  of 
possible  loading  fronts  art 


4>p  «  u 
♦S  "  77 


3K  +  4G 
o  o  ■ 


sin  (v  „  3f( 


.  -1  .1]  o, 
sin  —) 


(22) 


(23) 


In  the  present  problem  the  only  discontinuity  for  unloading¬ 
reloading  is  an  S  front  at 


4>s  -  n 


sin 


(24) 
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As  stated  previously,  only  superseismic  velocities  V 
will  be  considered,  i.e,,  the  value  of  V  will  be  sufficiently 
large  so  that  the  angles  $p  ,  ,  <£g  defined  by  Eqs.  (22,23,24) 

are  real  and  in  the  range  between  u/2  and  it.  In  addition,  due 
to  G  >  G  ,  the  inequality  <  <{..  must  hold  so  that  the 

0  0  O  o 

relative  locations  of  the  P,  S  and  S  fronts  are  shown  in 
Fig.  (4). 


Stress  Changes  at  $  =  <J>p 


The  first  stress  change  must  be  loading  because  the 

material  is  originally  -unstressed  and  unstrained.  Since 

<j)p  <  0 s  ’  the  first  stress  change  occurs  discontinuously  at 

4>  °  <(>p  .  The  relation  between  the  principal  st:  esses  0^ 

(normal  to  the  front)  and  ©2  (parallel  to  the  front)  behind 

the  P  front,  i.e.,  for  4>  >_  ,  is 

3K  -  2G 
o _ o 

°2  3K  +  4G  1 
o  o 

and  the  direction  of  the  principal  stress  0^  is 


(25) 


9  = 


(26) 


At  this  stage  the  value  0^  is  not  known  and  conditions 
for  its  determination  from  the  prescribed  surface  load  p  must 
be  developed. 


Stress  Changes  at  <J>  =  <J>g 

The  values  of  the  stresses  and  of  the  direction  9  for 
<$>  >  <}>p  must  remain  as  given  by  Eqs .  (25,26)  until  <{>  «=  4>g 


11 


where  a  second  fron..  of  discontinuity  may  occur.  The  stresses 
must  change  at_  this  S  front  in  such  a  way  that  the  material 
is  unloaded  and  possibly  subsequently  reloaded  as  indicated 
by  the  value  of'Jj  .  To  describe  the  S  front,  the  state  of 
stress  0 ^  ,  o 2  >  0  ahead  of  $  »  Is  resolved  into  components 
parallel  and  perpendicular  to  $  “  $g  and  the  resulting  com¬ 
ponents  of  stress  are  denoted  by  5^  ,  0^  and  x^  ,  Fig.  (5). 

Since  the  S  front  involves  only  shear  changes,  0^  and  oT  are 
constant  across  <}>  =  $  and  only  the  value  of  x.  can  change  to 
another  value  x0  .  However,  since  the  S  front  is  one  of 
unloading-reloading  the  value  of  for  x2  behind  the  front 
must  not  estceed  the  value  of  J2  for  T^  ahead.  This  implies 

-T i  <  T2  <  Tx  (27) 

Equation  (27)  has  a  geometric  interpretation  in  terms 
of  the  principal  directions  0^  ,  @2  on  either  side  of  the 
front.  Let  6  denote  the  smaller  of  the  two  angles  made  by 
the  normal  to  the  S  front  and  the  principal  directions  ahead  ** 
of  the  front.  (The  two  possible  cases  are  shown  in  Figs.  (6a) 
and  (6b).)  Equation  (27)  implies  that  the  appropriate 
principal  direction  behind  the  S  front  makes  an  angle  with 
the  normal  to  the  front  which  is  less  than  or  equal  to  6, 
i.e.,  it  lies  in  the  shaded  wedge  in  the  respective  Fig.  (6a) 
or  (6b).  If  this  shaded  wedge  includes  a  vertical  line,  Fig. 

(7),  a  solution  involving  only  the  P  front  and  S  front  is 
possible.  In  this  case  the  value  of  the  ratio  x2/x^  can 
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chosen  so  as  to  make  the  appropriate  principal  axis  behind 
the  S  front  vertical,  insuring  that  the  boundary  condition 
of  vanish'ng  shear  stress  is  satisfied.  The  value  of  can 
then  be  chosen  to  satisfy  tne  boundary  condition  0^  =  p. 


It  is  easy  to  derive  a  condition  indicating  when  a 
solution  of  the  type  just  described  is  possible.  Figure  (7) 
shows  the  P  and  S  fronts  and  the  shaded  areas  which  are  to 
contain  the  vertical.  If  Fig.  (6a)  applies,  the  relation 

-  <J>p  =  6  _<  ^  holds,  so  that  the  vertical  line  will  be  in¬ 
side  the  shaded  wedge  if  it  -  <_  6  =  <j>s  -  4>p  •  This  con¬ 

dition  may  be  written  in  the  form 

*  +  <S>P 

A  N. 


(28) 


If  Fig.  (6b)  applies  no  additional  condition  is  required 
because  the  wedge  always  contains  the  vertical  line.  Noting 
that  in  this  case  >  ~  and  5.  >  <>_  +  y  one  finds  that 
Eq.  (28)  is  automatically  satisfied.  It  is  therefore  the 
general  condition  for  the  existence  of  a  solution  involving 
only  a  P  front  and  an  S  front. 


If  Eq.  (28)  is  not  satisfied,  it  may  be  possible  to 
construct  solutions  by  including  at  $  E  <f>g  an  additional 
shear  front  at  which  loading,  occurs.  This  possibility, 
to  be  considered  next,  can  of  course  occur  only  if  Fig.  (6a) 
applies . 
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Stress  Changes  at  $  *■  (Jig 


Immediately  ahead  of  the  S  front,  Fig.  (8),  the  state  of 
stress  is  identical  to  that  immediately  behind  the  S  front. 
Since  the  S  front  inherently  involves  additional  loading,  the 
value  of  in  the  region  between  the  S  and  S  fronts  must  be 
equal  to  the  previous  maximum  value  of  J2  between  the  P  and  S 
f.ronts.  This  implies  that  t 2  “  at  the  S  front  as  well  as 

t2  >  T ^  at  the  S  front,  where  and  x2  are  t'le  components  of 
the  shear  stress  parallel  to  the  S  front  for  arid 

<J>  >  <J)  ,  respectively,  Fig.  (8).  Noting  the  definition 

b 

6  “  ^5  ■  4>p  •  the  first  reqjirement,  X2  “  -T^  ,  implies  that 
the  direction  02  of  the  principal  axis  behind  the  S  front  is 


5S  +  &  -  f  =  2$s  - 


TT 

2 


(29) 


This  is  equivalent  to  saying  that  the  principal  axis  coincides 
with  the  right  edge  of  che  shaded  wedge  in  Fig.  (6a).  The 
second  requirement,  x2  >  X^  ,  implies,  Fig.  (8),  that  the 
angle  between  the  normal  to  the  front  and  the  direction  of 
the  principal  stresses  ahead  of  and  behind  the  front, 

®2  “  ~  §)  and  "  4>g)  >  respectively,  must  satisfy  the 

inequality  82  -  (<J>S  -  <  n  -  <J) g  .  T.n  addition,  the 

direction  of  the  principal  stress  between  the  S  and  S  fronts 
must  be  at  least  as  steep  as  the  normal  to  the  S  Clout,  or 
®2  -  ~  §  *  Substitution  of  Eq.  (29)  into  the  last  two 

inequalities  gives 


♦s  + 


it  + 


1  *s 


2 


2 


(30) 
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IV  NONLINEAR  MATERIAL. 


A.  Material  Description 


The  fora  of  the  moduli  K  and  G  as  functions  of  the  in¬ 
variants  was  '.aken  from  Ref.  [2],  For  initial  loading, 

0  -  Go  -  I^l  -  Vl  +  Vl 


K  "  Ko  "  Klekk  +  Vkk 


(33) 

(34) 


and  for  unloading  and  reloading, 


G  +  G.  /J7  -  G,J.  +  G,jf 

O  1  l  l  1  J  X 


K  -  K,J, 
o  11 


(35) 

(36) 


The  magnitudes  and  signs  of  the  constants  are  to  be 
selected  such  that  the  stress-strain  diagrams  have  the 
curvatures  and  general  character  shown  in  Figs,  (la)  and 
(2a). 


B .  Method  of  Solution 

The  configurations  of  the  solutions  shown  in  Figs. 

(7,8)  for  the  bilinear  material  are  used  as  guides  in  the 
determination  of  solutions  in  the  nonlinear  case.  Beginning 
with  a  discontinuity  in  the  principal  stress  normal  to  the 
front,  $  =  ,  Eq.  (A-3),  it  can  be  expected  that  instead 

of  the  discontinuities  in  shear  in  Figs.  (7)  or  (8),  the 
nonlinear  material  will  produce  regions  in  which  shear  stresses 
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change  in  a  continuous  manner.  Because  of  the  change  in 
curvature  in  the  unloading-reloading  portion  of  the  diagrams 
Fig.  (2b),  there  is  also  a  possibility  of  a  noncharacteristic 
shock  in  shear  adjoining  the  continuous  region. 


Conditions  which  define  the  discontinuity  in  normal 
stress  at  <}>  =  are  discussed  in  Appendix  A.  Using  an 

p 

inverse  approach,  an  arbitrary  jump  in  the  principal  stress 
0^  normal  to  the  front  is  selected.  Integration  of  the  in¬ 
cremental  relation  at  the  front,  Eq.  (A-3)  yields  the  corre- 
P  P 

spending  values  o ^  and  at  the  front,  and  hence  the  starting 
values 


2  frr*  P  \ 

3  "  °2 


1  P 

s2  “  "  2  S1 


.  J, 


P  P 

‘i  +  202 


Fov  locations  <f>  >  <J>p  ,  Eqs .  (16)-(19)  govern  the 
continuous  portion  of  the  solution,  where  the  moduli  are 
given  by  Eqs.  (33)-(36).  Equations  (16)  may  be  satisfied 

iiit 

.by  regions  of  uniform  stress,  si  "  s2  “  Ji  "  ®  “  0 ,  or  by 

regions  of  stress  variation  when  Eq .  (19)  is  satisfied. 
Because  Eq.  (19)  implies  that  the  determinant  of  Eqs.  (16) 
vanishes,  only  three  of  the  four  Equations  (16)  are  thus 
independent.  However,  since  Eq .  (19)  must  remain  valid 


(37) 
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throughout  regions  of  stress  change,  a  fourth  independent 

i  •  i  t 

equation  in  Sj^  ,  s,,  ,  ^  and  ^  -  s2>0  is  obtained  by 
differentiation  of  Eq .  (19)  with  respect  to  $.  In  the 
present  problem  the  resulting  relation  is  of  interest  only 
for  the  case  of  changes  in  the  deviatoric  stresses,  when 

t 

■  0.  For  this  case  the  fourth  equation  becomes 


(2sl  + 


Vs!  + 


(2s2  +  s^ 


2pV 


sin  2<j> 


(38) 


*) 


Substitution  of  Eqs .  (21)  and  (38)  into  Eqs.  (16)  yields 
the  following  expressions  for  the  four  unknown  derivatives 

'2  I/J7I 

sin  2<f> 


2pV 


G1^S1  ”  s2 


(S1  “  s2^  Can 


(?V> 

(40) 

(41) 

(42) 


*) 


Equations  (39)— (42)  permit  the  numerical  determination  of 
the  values  of  stresses  in  regions  of  continuous  loading  or 
unloading-reloading  by  quadratures,  provided  the  values  on 
one  boundary  of  the  region  are  known. 


The  possibility  of  a  noncharacteristic  shock  in  shear 
associated  with  the  solution  must  also  be  explored.  From 
jKg  case  of  uni 3x1  ul  pressure  VwVcs  i. s  vcl.1  knovn 

in  a  material  which  continuously  hardens  as  the  stress  level 


The  value  G.^  shown  applies  for  the  case  of  loading. 

For  unloading  or  reloading  it  is  to  be  replaced  by  . 
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is  changed  from  level  "A '*  to  a  level  "B",  discontinuities 
with  a  jump  in  stress  from  "A"  to  "B"  can  occur.  The  stress- 
strain  diagram  for  shear,  Fig.  (9),  indicates  that  this  is 
the  situation  in  the  reloading  portion,  between  points  3  and 
4  of  the  diagram.  The  complete  unloading-reloading  diagram 
has  a  point  of  contraf lexure  at  3.  In  such  a  case  it  is  also 
known  from  prior  experience  that  there  will  be  a  continuous 
stress  change,  immediately  followed  by  a  discontinuity. 

The  discontinuity  and  its  velocity  can  be  represented 
graphically  using  the  stress-strain  diagram.  Fig.  (9).  Stress 
and  strain  may  change  discontinuously  from  a  value  corresponding 
to  a  point  2a,  between  2  and  3,  to  a  point  3a,  between  3  and  4‘. 
Stability  of  the  shock  requires  that  a  line  3a  to  2a  is  tangent 
to  the  T-y  diagram  at  point  2a.  The  slope  of  the  line  defines 
the  shock  velocity.  As  the  diagram  applies  only  as  far  as 
point  4  the  stress  level  cannot  exceed  the  value  at  point  4. 

If  the  discontinuity  present  in  a  solution  does  not  reach 
a  stress  level  corresponding  to  point  4  where  the  value  of  J2 
is  equal  to  the  previous  peak,  the  solution  can  not  have  a 
region  of  further  loading  in  shear.  If  the  discontinuity  is 
such  that  a  value  ,  equal  to  the  previous  peak,  is  reached 
a  region  of  further  loading  is  possible. 

C .  Determination  and  Description  of  the  Solution . 

Figure  (10)  shows  the  configuration  of  the  solution  which 
will  be  constructed.  It  is  a  generalisation  of  the  one  in 


Fig.  (8)  for  the  bilinear  material  in  which  two  discontinuities 
in  shear  occur,  one  in  unloading-reloading,  one  in  further 
loading . 

At  the  invariant  J,  changes  discontinuously  from  J„  ■  0 

to  a  value  •  The  direction  of  the  principal  stress  is 

P  TT  P 

defined  by  the  angle  0  a  ^p  ~  J  »  while  the  value  of  02  corre- 

p 

sponding  to  the  chosen  value  can  be  computed  from  the  re¬ 
lations  given  in  Appendix  A.  Behind  the  front,  $><{>  ,  there 

is  a  region  of  uniform  stress  which  terminates  at  the  root  <J>^ 
of,  Eq.  (21)  using  the  value  of  G  for  unloading.  At  this 
location  a  region  of  continuous  unloading  of  deviatoric  and 
shear  stresses  begins.  Forward  integration  of  Eqs.  (39)-(42)-, 
using  the  initial  values,  Eq.  (37),  determines  the  stress 
field  numerically,  point  by  point. 

The  integration  is  stopped  at  a  trial  location  <t>2  > 
which  is  selected  so  that  the  value  of  0,  defining  the 
direction  of  the  principal  stress,  satisfies  the  inequality 
0  <  <P  -  ^  •  (This  ensures  that  a  stable  shock  in  shear  in  the 
location  <j>2  Is  possible.)  As  next  step,  the  location  <J>g  of  a 
shock  in  shear  is  determined  in  which  the  shear  stress  changes 
from  the  value  found  at  $2  to  the  value  which  makes  J2  equal 
to  the  previous  maximum  of  J2  .  The  point  <{>2  and  the  shock 
described  are  accepted  if,  and  only  if,  <}>g  =  <j>2  ,  otherwise 
another  point  <J>2  must  be  selected,  and  the  process  repeated. 

For  <f>  >  $2  =  $5  there  is  again  a  region  in  which  the 
principal  stresses  remain  constant.  This  region  can  be 


—  20  -- 


terminated  at  a  root  $  **  $3  of  Eq .  (21)  using  G  for  loading, 
in  which  location  a  new  region  of  continuous  stress  change 
begins.  Forward  integration  of  Eqs.  (39)-(42)  proceeds  until, 
at  (j>  ■>  ,  the  principal  stress  0^  becomes  perpendicular  to 

the  surface,  i.e.,  0  *  -g  •  From  the  point  <{>  **  <j>i  to  the 
surface,  <}>  »  it  ,  the  stresses  are  again  uniform.  The  value 
of  the  principal  stress  obtained  in  the  location  § ^  is  the 
value  of  the  applied  surface  pressure  which  corresponds  to  the 

P 

originally  assumed  value  0^  . 

There  is  also  the  possibility  of  a  solution  which  does 
not  contain  a  region  of  change  associated  with  further  loading 
in  stress  between  4> ^  and  .  This  solution  corresponds  to 
the  one  in  Fig.  (7)  in  the  bilinear  case.  In  this  solution  a 
shock  can  be  selected  such  that  the  principal  stress  for 
(j>  >  <j>2  is  vertical. 


(43) 


(44) 


D .  Numerical  Results 

The  moduli  chosen  for  the  first  two.mumerical  examples 

are 

G  =  7.178  -  20.0  | /I7| 

K  »  11.965  -  150.0  ekk 

for  loading,  and 

G  *  4.0  +  40.0  |/J^| 


1 


for  unloading  and  reloading. 
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Because  the  modulus  K  affects  tha  solution  only  during 
loading  at  4>  -  <J)p  no  information  for  unloading  is  required. 

It  may  be  seen  from  Eq .  (39)  that  solutions  in  regions  of 
continuous  stress  changes  (in  shear)  are  influenced  only  by 
coefficients  and  .  Without  actual  loss  of  generality 
G^  ,  Gj  ,  Gj  ,  G^  were  conveniently  assumed  to  be  zero. 

Figures  (11)  and~(12)  show  numerical  results  for  the 

moduli  given  above  for  two  different  ratios  of  surface  pressure 

speed,  V,  to  the  speed  cp  .  It  may  be  noted  that  the  regions 

of  continuous  stress  change  are  quite  small  in  angular  extent. 

This  reflects  the  relatively  small  amount  of  curvature  of  the 

stress-strain  diagrams  for  the  moduli  used.  There  is  a 

relatively  small  change  in  0^  from  the  value  at  <{>  “  <f>?  to  the 

surface  value  at  $  ■  ir.  However,  if  stress  components  in  the 

horizontal  and  vertical  directions,  a  and  <J  respectively, 

x  y 

are  considered  [see  Figs.  (13)  and  (14)]  the  changes  In  the 
stress  field  are  more  significant  and  visible. 

Figure  (15)  shows  numerical  results  for  the  modulus 

G  «  4.0  +  10.0  j/j^j  <*5) 

for  unloading  and  reloading.  The  initial  loading  moduli  were 
chosen  to  be  the  same  as  Eqs.  (43).  In  this  case  only  one 
region  of  continuous  stress  change  is  present.  The  stress 
changes  introduced  by  the  shock  in  shear  are  sufficient  to 
satisfy  the  boundary  conditions,  and  no  region  of  further 
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loading  is  -equired.  This  alternate  solution,  corresponds 
to  the  one  shown  in  Fig.  (7)  for  a  bilinear  material.  It 
occurs  if  the  parameters  are  such  that  a  shock  in  shear  can 
be  selected  so  that  the  principal  stress  becomes  vertical, 
0  "  —  ,  immediately  behind  the  front.  In  this  case  the  front 
is  represented  by  the  line  2a-3a  in  Fig.  (9). 
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APPENDIX  A 


Consider  a  plane  moving  front  at  which  discontinuities  in 
stresses  and  velocities  occur,  and  let  £,ri,C  be  a  Cartesian 
coordinate  system,  Fig.  (A-l),  where  5  is  normal,  and  r) ,  5 
tangential  to  the  front.  Conditions  for  discontinuities  of 
stress  and  velocity  normal  as  well  as  tangential  to  the  front 
will  be  studied,  but  the  two  cases  must  be  treated  separately. 


1 .  The  P-front. 

A  front  where  discontinuities  of  stress  and  velocity 
occur  in  the  direction  5  of  propagation  of  the  front  will  be 
referred  to  as  a  P-front.  Let  Ao^  ,  Av^  be  the  discontinuous 
changes  of  normal  stress  and  velocity,  respectively.  Conser¬ 
vation  of  momentum  at  the  front  requires 

A <*£  “  pCpAv^  (A-l) 


where  cp  is  the  velocity  of  propagation  of  the  front,  while 

\ 

<\er  «■  —  A v» 

?  cp  C 


Av„  =  Av_  -  0 

n  C 


Ao. 


Ao^  i  0 


(A-2) 


J 


Equation  (A-l)  and  the  first  of  Eqs .  (A-2)  give  the  relation 
o  Ao - 


pAc, 


The  inclination  <*'p  of  such  a  front  in  the  problem 


considered  in  the  body  of  the  report  is  thus 


.  -1  rl 
tt  -  sin  [y 


Ao, 


pAe 


( A—  3) 


C 
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Consider  in  detail  the  case  where  the  region  ahead  of  the 
front  is  stressless  and  at  rest.  The  normal  to  the  front 
is  then  a  principal  direction,  i.e.,  *  a2  “  ar\  * 

.  To  find  for  any  selected  value  Ao^  the  jumps  in 
the  other  quantities  Eqs.  (1)  and  (2)  are  combined  with 
Eqs.  (A-l)  and  (A-2)  to  obtain  the  "one  dimensional"  incre¬ 
mental  stress-strain  relation 

Oi  -  (K  +  ~  O 


Using  the  values  of  K  and  G  given  by  Eqs.  (33)  and  (34)  the 
integral  of  d 2  can  be  found  in  closed  form.  First  Eq.  (1) 
is  integrated  to  express  as  a  function  of  •  'The 

result  is  substituted  into  Eqs.  (5),  (6)  and  (33),  then 
into  Eq.  (A-4) ,  with  the  result 


where 


do  „ 

d^'+  a02 


l  Cj  ei 
J-l  J 


J-l 


f  /3  Of 


K  -  -5-  G 
o  3  o 


C2  -  -Kx  +  aKo  +  2G2Ko 


K2  -  j  aKj  -  G2K1  -  6G,K 


2 
3  o 


C4  "  1  aK2  +  3  C2K2  +  6G3KoK1 


Cc  -  -4K  K,G, 
5  o  2  3 


C,  -  2G3K1K2 
C7  -  "  !  G3K2 


3  2 

i  g3  x 


(A-4) 


(A-5) 


(A-6) 


» 
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The  solution  of  the  differential  equation  (A-5) ,  for  the 
initial  condition  a 2  =  0  when  £  ^  «=  0,  is 


7  7 


>2  -  -2  l  l 

m=l  J>=m  a 


(n-1) ! 


C . (1-6.  e 
3  lm 


•aCl)  er1  T  3K  £ 
1  0 


-1.5  K1e^  +  K  ( A-7 ) 


where  6,  “1  for  m=*l,  6,  =0  for  m  t7'  1. 

lm  ’  lm 


A  Newton  Raphson  technique  was  used  in  the  examples  to 


find  the  value  of  for  the  selected  values  of  0 


a  may 


then  be  found  from  the  integral  of  Eq .  (1)  giving  the  complete 

stress  system  at  the  front. 

2.  The  S-front. 

Discontinuous  changes  at  the  S-front  occur  only  in  the 
velocity  v^  parallel  to  the  front,  and  in  the  corresponding 
stress  a.  **  T.  Conservation  of  momentum  at  the  front 
requires 


Ax  «  pcsAv^ 


(A-8) 


where  Cg  is  the  velocity  of  propagation  of  the  front,  while 


ay 

Av 

Ao 


- —  Av 
2c ,  n 

d 

t  J  av?  '  0 


AO. 


Ao?  =  0 


(A-9  ) 


Equation  (A-8)  and  the  first  of  Eqs.  (A-9)  give  the  relation 
2  At 

CS  "*  TpAy  ’  T^e  inclination  g  of  such  a  front  in  the  problem 
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considered  in  Che  body  of  the  report  is  thus 


4>s  *=  IT 


s  in 


-1 


(A-10) 


To  determine  the  ratio  Ai/Ay  consider  the  case  where  the 
front  propagates  into  a  region  which  has  an  existing  stress 
field  with  principal  directions  not  perpendicular  to  the 
front.  The  incremental  stress-strain  relation  is 


T  •=  2G  Y 


(A-ll) 


where  0  depends  on  the  invariant 

J2  “  sc  +  s£sn  +  sn  +  t2  < a— 12) 

As  stated  in  the  body  of  the  report,  only  discontinuities  in 
shear  for  unloading  -  reloading  were  required.  For  this 
situation  Eq.  (A-ll)  was  integrated  numerically  using  Eq.  (35) 


for  G. 
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APPENDIX  B 


Bilinear  Fluid. 

The  simplest  case  of  the  new  material  model  considered 
in  this  report  is  that  of  a  bilinear  fluid  with  a  pressure- 
volume  relation  similar  to  the  one  shown  in  Fig.  (lb)  for 
the  solid.  If  a  half-space  of  such  a  material  is  subjected 
to  a  supersonic  progressing  step  load,  the  solution  is  trivial. 
It  consists  of  a  shock  front  (in  which  the  lower  modulus 
for  loading  applies)  follr  ed  by  a  uniform  stress  field.  The 
higher  value  of  the  modulus,  a  ■'s  not  affect  the  solution 
at  all.  The  simpler  nature  of  the  bilinear  fluid,  as  com¬ 
pared  to  the  solid,  makes  it  possible  to  give  a  closed  form 
solution  for  the  alternative  problem  of  a  half-space,  Fig. 

(B-l) ,  subject  to  a  progressing  pressure  pulse  which  rises 
suddenly  to  a  value  ,  but  subsequently  decreases  linearly. 
This  problem  involves  both  loading  and  unloading  in  the  fluid, 
so  that  the  nature  of  the  material  affects  the  result. 

The  basic  equations  describing  an  inviscid  fluid  may  be 
written  in  nonaimensional  form  by  introducing  the  following 
units  of  length,  veloci'.  ,  and  pressure,  respectively, 

p 

_ o  _  discontinuity  in  surface  pressure 
l  gradient  of  surface  pressure 

V  -  velocity  of  surface  load 

»oy2 

where  p  is  the  nominal  density  of  the  fluid.  The  linearized 

o 

nondime n sion a  1  equations  of  conservation  of  mass  and  moment  urn 


length : 

veloci ty : 
press  ure : 


* 
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are 


«  0 

3x  dt 

If  +  fe  - 0  >  <B'1) 

H+  (-T?>  <17  -  If5  ■  0 

p  V  / 

O 

where  p  is  the  pressure  in  the  fluid  while  u  and  v  are  the 

horizontal  and  vertical  components  of  velocity,  respectively. 

In  this  form  the  equations  depend  solely  on  the  nondimens ional 
£ 

parameter  - r  .  This  parameter  varies  because  K  assumes  the 

poV 

value  or  depending  upon  whether  the  fluid  is  being 
loaded  or  unloaded  -  reloaded,  respectively. 


In  the  steady-state  problem  considered  here,  p,  u  and  v 
ar-  necessarily  functions  of  y  and  the  combination  5  ■*  x-t. 
After  elimination  of  u,  Eqs .  (B-l)  can  therefore  be  reduced  to 


jjp  3v  r 

3?  '  95  “  C 

2 2  l£  _ 

9C  3y 


(B-2) 


where  Z  is  the  nondimensional  parameter, 


poV 


(B-3) 


’which  may  assume  the  alternative  values  Z^  ,  Zy  for  loading 
and  unloading  -  reloading,  respectively. 


Consider  the  situation  when  V  is  supersonic  with  respect 
to  the  sound  speeds  for  loading,  [/  ~  ,  and  unloading, 


A3 


If  the  fluid  were  linear  with  K  =  K  ,  instead 

L 

of  bilinear,  one  finds  easily  that  the  solution  consists  of  a 

shock  front  inclined  at  the  angle  4>p  •*  it  —  sin 

followed  by  gradually  (as  a  matter  of  fact  linearly)  decreasing 

pressure.  It  will  be  assumed,  and  subsequently  verified,  that 
the  situation  for  the  bixinear  fluid  is  similar,  i.e.,  there 
is  loading  at  the  front  only,  and  unloading  everywhere  else. 

As  the  entire  region  between  the  surface  and  the  ?-front, 

Fig.  (B-2) ,  is  assumed  to  be  an  unloading  region,  the  problem 

requires  the  solution  to  the  hyperbolic  differential  equations 

(B-2)  for  Z  =  Zy  in  this  region.  This  is  now  a  linear  hyper-' 

bolic  problem.  As  the  boundary  conditions  for  this  problem, 

one  imposes  at  all  points  Q  on  the  loading  front,  i.e., 

C  *  -Z  y  ,  the  momentum  requirement 
x  "  x 

vQ  -  ZLpQ  (B-A) 

while  at  any  surface  point  T  with  the  coordinate  £  ,  the 

surface  pressure  is  specified, 

PT  “  po(1  +  ?T)  (B-5) 

The  linear  form  of  the  boundary  conditions,  Eqs .  (B-A) 
and  (3-5),  in  this  linear  hyperbolic  problem  implies  that  the 
pressure  p (£ ,y )  is  of  the  linear  form 

p(5»y)  =  P  0  C 1  +  5  +  ky)  (B-6) 

where  the  coefficient  of  £  is  unity  in  order  to  satisfy 
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Eq.  (B-5),  while  the  constant  k  is  so  far  undetermined . 
Substitution  of  Eq.  (B-6)  into  Eqs.  (B-2)  leads  to 


■v(i;,y)  ■=  p0(vq  +  k£  +  Zy y ) 


(B-7) 


where  vq  is  a  constant  of  integration.  The  values  of  k  and 
vq  are  determined  by  substitution  of  Eqs.  (B-6)  and  (B-7) 
into  Eq.  (B-4) , 


k 


v 

o 


Z 


L 


(B-8) 


•In  order  to  verify  the  premise  that  the  entire  wedge 
between  the  P-front  and  the  surface  is  continuously  unloading 
consider  the  pressure,  obtained  from  Eq.  (B-6),  p  =  pQ(l  +  x  “ 
t  +  ky).  Since  is  obviously  negative,  the  validity  of  the 
solution  is  confirmed.  Equations  (B-6),  (B-7),  (B-8)  and 
Fig.  (B-2)  thus  describe  the  solution  to  the  problem. 


The  solution  presented  above,  while  intended  for  com¬ 
pletely  supersonic  velocities  V  of  the  surface  pressure, 
appears  to  remain  valid  when  the  load  moves  subsonlcally  with 


respect  to  the  higher  of  the  sound  speeds 


fir 

•ip' 


provided 


that  V  is  still  supersonic  with  respect  to  the  sound  speed  in 

i  [si 

loading,  /  — —  .  It  is  of  general  interest  to  show  that  in 
po 

this  case,  where  the  situation  is  not  completely  supersonic, 


t^e  above  solution,  while  valid,  is  not  unique.  Alternate 
solutions  are  obtained  by  adding  terms  which  satisfy  homogeneous 
boundary  conditions  on  the  surface  and  on  the  P-front. 
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For  the  case  now  being  considered,  Eq .  (B-3)  gives  a 
2  2 

negative  value  for  Z  while  Z.  remains  positive.  Introduce 

U  L 

the  real  and  positive  quantity  Y, 


2  2 
Y  -  -Z„ 


Y  >  0 


(B-9) 


and  a  new  set  of  independent  variables, 


r  =  ]/  ^  +  y2 
6  =  tan  1  (-  ^ 


) 


(B-10) 


It  may  be  verified  easily  that 


1  2<f> 

pa— 

Y  o  y 

v 

v  =  -y  n 


(B-ll) 


with 


<J>  =  rA  cos  A0 


( B-12 ) 


is  a  solution  of  Eqs.  ( B— 2 )  which  satisfies  the  homogeneous 
boundary  condition  p  =  0  on  the  surface  0=0  for  all  A. 

The  constant  A  can  be  used  to  satisfy  the  r emaining  conditions, 

Eq.  (B-4).  If  singularities  in  pressure  and  velocity  are 
eliminated  as  inappropriate,  one  obtains  an  infinite  number 
of  roots  A  =  Ar  ,  namely, 

A  -  - ;  n  =  1,  2,  ...  (B-l  3) 

4l 

where  the  principal  value  of  tan  ^  is  used.  Therefore,  one 


* 
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may  add  to  the  previous  solution,  Eqs.  (B-6) ,  (B-7),  (B-8)  and 
(B-9),  a  term  of  the  form  of  Eqs.  (3-11)  where 

CO 

<j>  «  T  A  r^n^cos  X  6  (B-14) 

,  n  n 

n=l 

Tne  coefficients  are  arbitrary,  except  for  the  restriction 
that  ^  <  0  everywhere.  Steady-state  solutions  are  thus  not 
unique  when  the  surface  load  moves  subsonically  with  respect 
to  the  speed  of  sound  which  applies  for  unloading  -  reloading. 
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